Interplay of Rayleigh and Peierls Instabilities in Metallic Nanowires 
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A quantum-mechanical stability analysis of metallic nanowires within the free-electron model 
is presented. The stability is determined by an interplay of electron-shell effects, the Rayleigh 
instability due to surface tension, and the Peierls instability. Although the latter effect limits the 
maximum length also for wires with "magic radii", it is found that nanowires in the micrometer 
range can be stable at room temperature. 
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Metal nanowires represent the ultimate limit of con- 
ductors down to a single atom in thickness. In recent 
years, experimental research on metal nanowires has bur- 
geoned and various techniques have been employed to 
produce them for many different metals [1, 2, 3, 4, 5]. 
Remarkably, these nanowires are quite stable, although 
one might expect them to break up because of surface 
tension. Such an instability arises from classical contin- 
uum mechanics and is know as the Rayleigh instability 
[6]. Experimentally the stability was found to depend 
sensitively on the geometry of the wire. Yanson et al. 
[2, 3] have measured conductance histograms for alkali 
metal nanowires. At low temperature only wires with 
certain "magic radii" were stable. In another notable 
experiment [4, 5], Kondo et al. have produced long gold 
nanowires that were almost perfectly cylindrical in shape. 
Electron microscope images show lengths between 3 and 
15 nm and diameters from 0.6 to 1.5 nm. A histogram 
was built which again reflects stability for certain diam- 
eters only. 

Many features of metallic nanowires can be explained 
in terms of a model of free electrons confined to the ge- 
ometry of the wire [7]. At least for alkali metals, exper- 
iments on conductance quantization and current noise 
are in good agreement with theoretical predictions based 
on the Jellium model [8, 9, 10, 11, 12]. Furthermore, 
Kassubek et al. [13] have shown that a semi-classical 
calculation of the density of states of nanowires leads 
to stable magic radii, like those observed in experiments 
[2, 3]. However, a puzzling conclusion of the semiclassi- 
cal theory is that nanowires with a magic radius remain 
stable whatever their length. 

There is substantial interest in the stability of 
nanowires, since they may serve as conductors in fu- 
ture nanocircuits. In particular, one would like to know 
whether they remain stable up to a length useful for ap- 
plications. We therefore reexamine this problem using a 
fully quantum-mechanical analysis. It turns out that the 
suppression of the Rayleigh instability by electron shell 
effects, which is captured by the semiclassical theory, is 
now supplemented by an interplay between the Rayleigh 
and a novel Peierls-type instability missed in previous 
work. In fact, this latter quantum mechanical instability 



limits the maximal length of stable nanowires but still al- 
lows for wires in the micrometer range sufficient for most 
applications. 

We treat the metallic nanowire as an open system con- 
nected to macroscopic metallic electrodes at each end. It 
is most naturally described within a scattering matrix 
approach [9]. The relevant thermodynamic potential is 
the grand canonical potential which is related to the 
density of states D{E) by 
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where fcs is the Boltzmann constant, T is the temper- 
ature and /i is the chemical potential specified by the 
macroscopic electrodes. The density of states can be cal- 
culated from the scattering matrix S{E) by 



(2) 



A factor of 2 for spin degeneracy has been included. In or- 
der to examine the stability of a nanowire, we determine 
the energy change (SSI caused by a small deformation of 
the wire geometry. A nanowire is only stable if all small 
deformations result in an increase of 57. 

The Jellium model considers free electrons that are 
kept in the wire by a confining potential chosen as a hard- 
wall boundary. The corresponding Schrodinger equation 
can be solved by a separation ansatz 
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in cylindrical coordinates r, and z. This is motivated 
by the fact, that at a given position z along the wire we 
can find a complete set of orthogonal transverse eigen- 
functions Xfj-i"!^, 'P', z) with eigenenergies E^{z). 

Inserting this ansatz into the Schrodinger equation, we 
obtain a system of coupled differential equations for the 
coefficient functions Cf_,{z) [14] 
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where a prime indicates differentiation with respect to z. 
For convenience, we use units with fi^ /2me = 1, where 
rrif, is the electron mass. The coupling matrix reads 
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In the following we want to examine a cylindrical wire 
and its sole classically unstable deformation - an axisym- 
metric one. In this case the geometry can be character- 
ized by the radius function R{z) and the coupling matrix 
can be written as A^^{z) = Af_i,^ ■R'{z)/R{z), where Afj.i, 
is now independent of z. Eq. (4) simplifies to read 
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In the adiabatic approximation, which is justified if the 
radius varies slowly with z, the right hand side of Eq. (6) 
is neglected and the functions Cf^^z) can be determined 
from a WKB approximation [9]. This procedure is not 
valid for small wavelength modulations of the wire ge- 
ometry that will be shown to limit the maximal length. 
Hence, we have to solve Eq. (6) systematically. 

To determine the stability of a cylindrical wire with 
given radius Rq and length L, we study small perturba- 
tions of its geometry given by 
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The Fourier coefficients bq are chosen such that the to- 
tal volume of the wire is unchanged by the deformation. 
Other physically reasonable constraints [15] are possible, 
but lead to similar results. We calculate the scatter- 
ing matrix of the deformed nanowire connected to two 
cylindrical leads by matching solutions of the Schrodinger 
equation (6) at the boundaries. The density of states (2) 
and the grand canonical potential (1) are then derived as 
an expansion in the deformation parameter s. 
The grand canonical potential is found to read 
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The leading order change il^^' due to the deformation is 
of second order in the Fourier coefhcients bq and can be 
written as 
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where the terms of order 1 / L include nondiagonal contri- 
butions which can be neglected if the wire is long enough. 
Hence the function a{q, Rq, L,T) essentially determines 




FIG. 1: a) Stability diagram for a cylindrical wire of length 
Lkp = 1000 at zero temperature. Shaded areas are unstable, 
b) unstable radii as a function of the wire length 



the stability of the nanowire and is therefore called sta- 
bility coefficient. If a is positive for any value of the 
perturbation wave vector q, all small deformations of the 
cylindrical geometry lead to an increase of Q and are 
therefore suppressed. 

First we want to discuss the limit of zero tempera- 
ture. Fig. 1(a) shows the stability diagram for a wire of 
length Lkp — 1000 as a function of radius Rq and per- 
turbation wave vector g at T = 0, where kp is the Fermi 
wave vector. Shaded areas indicate a negative value of 
a and therefore instability. For comparison, the dashed 
line shows the stability criterion for the Rayleigh insta- 
bility {qRo = 1), and the dotted lines show the criterion 
of the Peierls instability 
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where kp "' is the Fermi wave vector of band i^. 

A sequence of stable and unstable radii can be ex- 
tracted from such a stability diagram. This reflects the 
existence of so called magic stable radii, like they have 
been observed in experiments. Fig. 1(b) shows the extent 
of the regions of instability as a function of wire length, 
discussed below. 

In order to understand this result it is useful to distin- 
guish two contributions to a (Fig. 2): 
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The first term describes quantum effects due to the chan- 
nel structure of electronic transport and the second one 
includes effects of surface tension and curvature energy. 

For a given wire radius, a'"^' as a function of q shows a 
series of distinct negative minima. With increasing wire 
length, these minima get sharper and deeper, resulting in 
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FIG. 3: Dispersion relation of a channel v (solid line) due 

I to a perturbation of the wire with a wave vector q = 2fe^\ 

The inset shows the whole curve, while the main figure en- 
larges the vicinity of fe = fel^' . The dashed line represents the 



FIG. 2: The stability coefficients a. (soUd Une), a^"'^^ (dashed 
Une) and a^^"''^) (dotted line) as a function of q for R^kp = 
5.75 and Lkp = 1000. Inset: a^''^'> as a function of radius 
at g = O.lfeF and temperature T = (dotted line) and T = 
0.02Tf (solid line), Tp being the Fermi temperature. 



a divergence for an infinitely long wire. These instabili- 
ties arise from an effect similar to the Peierls instability 
[16]. Nanowires are not purely one dimensional, but due 
to quantization perpendicular to the wire axis, electron 
transport is divided into distinct channels. Each chan- 
nel has a quadratic dispersion relation and starts to con- 
tribute at a certain threshold energy, i.e., the cigcnencrgy 
of the corresponding transverse mode. This results in 
a sequence of quasi one dimensional systems with Fermi 
wave vectors k^p \ as defined in Eq. (10). A perturbation 
of the translational symmetry of the system results in the 
opening of energy gaps in the electronic band structure. 
If a gap opens at the Fermi surface, the total energy of 
the system is lowered. 

The detailed calculation here is different from the stan- 
dard theory of the Peierls instability of a perturbed linear 
chain with periodic boundary conditions [16]. There, the 
energy change increases linearly with the length L and 
a therefore is independent of L. Since the nanowire is 
part of a much larger system including the leads, the 
longitudinal wave vectors fc'^'^-' in the subbands are not 
restricted to multiples of ^ , like it is the case for an iso- 
lated system. Instead, the entire system formed by the 
leads and the wire determines the boundary conditions 
for the k^^\ Hence, an analogous Peierls problem arises if 
we consider a linear chain of length La that is perturbed 
by a potential ^ V cos(2A;f-z) acting within an interval of 
length L only, where L ^ Lq. In this case, the change of 
the eigcnencrgies of states near ifc^? can be calculated by 
perturbation theory and the resulting energy dispersion 
relation is sketched in Fig. 3. The dispersion relation is 
smooth. Only in the limit L ^ oo we recover the well 
known Peierls results with a jump at kp\ 



unperturbed parabolic dispersion. 



In the case of a single conducting channel, labeled v, 
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and has a minimum for q = 2k . For several open chan- 
nels, a'"'*' includes terms analogous to Eq. (12) for each 
channel, as well as similar but more complicated terms 
that take into account coupling between any two channels 
having the same azimuthal symmetry. This leads to addi- 
tional minima for q ^ kP+kP, that are less pronounced 
than those for q = 2kP . Away from the Peierls- type in- 
stabilities, Q*"*'' is positive and the nanowire is stabilized 
by electron shell effects. 

The second term in Eq. (11) also tends to stabilize the 
wire. It is given by 



q\ (13) 



with sums running over all open channels. A proportion- 
ality to q^ is expected for the change of the surface and 
curvature energy of the wire due to a deformation with 
wave vector g > 0. Using the explicit expressions for the 
coupling matrix elements one can show that a^""^'^ (g) is 
always positive. Note that the negative surface contri- 
bution a.t q = that leads to the well-known Rayleigh 
instability is contained in a.'^'^^\ 

Now a coherent description of the effects determining 
the stability of a nanowire arises. Surface tension and 
curvature energy tend to stabilize the nanowire for de- 
formations of small wavelength since the cost in energy 
increases oc q^ . On the other hand, the classically ex- 
pected Rayleigh instability for long wavelength perturba- 
tions is suppressed by electron shell effects that depend 
strongly on the wire radius. The channel structure of the 
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FIG. 4: (a) Stability diagram for a cylindrical wire of length 
Lkp = 1000 at three different temperatures. The shaded 
area indicates instability at T = O.OSTf, the solid and the 
dotted lines show the contours of the unstable regions for T = 
O.OITf and T = 0.005Tf respectively, Tf being the Fermi 
temperature, (b) unstable radii as a function of temperature 



tions to experiments. 

Recently experiments found evidence of nanowires 
made out of dimerized chains of gold atoms [18]. This 
is not in contrast to our results, as in this experiment 
surrounding molecular hydrogen had a stabilizing effect 
on the wire. Also the Jellium model used for our calcu- 
lations is at its limit of validity for a chain of atoms. 

In conclusion, we have presented a fully quantum me- 
chanical stability analysis of cylindrical nanowires, using 
the Jellium model. The stability of a nanowire is deter- 
mined by the interplay of surface tension, shell effects and 
the Peierls instability. Our calculated stability diagrams 
prove the existence of certain stable magic radii like they 
have been found in experiments and give an explanation 
for the observed cylindrical geometry of these nanowires. 

We thank C. A. Stafford, F. Kassubek and J. Biirki 
for valuable discussions. This research has been sup- 
ported by the DFG through SFB 276 and the EU Net- 
work DIENOW. 



electronic transport results in several Peierls-type insta- 
bilities for perturbations of wave vector q = 2fc^^ that 
become more pronounced with increasing wire length. 
The competition of these instabilities with the stabilizing 
effects of surface tension and curvature energy limits the 
maximal length of stable nanowires (Fig. 1). Neverthe- 
less, stable wire geometries persist up to lengths greater 
than Lkp = 10"* (i.e. L ~ Ifim for Na and Au) which is 
encouraging for the technological prospects of nanowires. 

Now we turn to the effect of temperature. Increas- 
ing temperature results in a broadening and a flattening 
of the minima describing the Peierls instabilities (Inset 
Fig. 2). The stability diagram of a cylindrical wire of 
length Lkp = 1000 is shown in Fig. 4(a) for three differ- 
ent temperatures. The Peierls instabilities disappear at 
higher temperature and the stability diagram approaches 
that of the classical Rayleigh instability, which is shown 
by the dashed line. Nevertheless islands of stability per- 
sist up to very high temperatures. Fig. 4(b) shows the 
extent of the unstable regions as a function of temper- 
ature and the increase of unstable configurations with 
growing temperature can be seen. The variation in the 
maximum temperature up to which the wires remain sta- 
ble reflects the supershell structure observed in Ref. [3]. 

In our calculations the free electrons are kept in the 
wire by an infinitely high potential well (" hard- wall" ) . 
Since the electrons in a real metal are not so strongly 
confined, a more physical boundary condition would be a 
potential well of finite height depending on the electron 
work function. The main influence of this "soft- wall" 
boundary condition is to increase the effective radius of 
the cylindrical wire relative to the "hard- wall" case [17]. 
This has to be kept in mind when comparing our predic- 
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